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' (-h '■ Abstract. For a general- form polarization biphoton qutrit, physically corresponding 

| to a pair of arbitrarily polarized photons in a single frequency and wavevector mode, 

we explicitly find polarization Schmidt modes. A simple method is suggested for 
factorizing the state vector and the explicit expressions for the factorizing photon 
creation operators are found. The degrees of entanglement and polarization of a qutrit 
are shown to depend directly on the commutation features of the factorizing operators. 
Clear graphic representations for the Stokes vectors of the qutrit state as a whole, its 
Schmidt modes, and factorizing single-photon creation operators are given, based on 
the Poincare sphere. 
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1. Introduction 

Biphoton qutrits are the simplest two-photon formations finding applications in the 
science of quantum information. In polarization biphoton qutrits two photons are 
assumed to have only one degree of freedom, polarization, with identical spatial and 
spectral features. Such states are usually produced via collinear spontaneous parametric 
down-conversion (SPDC) degenerate in the frequencies of the photons. The general form 
of a biphoton qutrit state vector is given by the superposition 

|*) = d|2tf) + C 2 \1 H , ly) + C 3 |2y), (1.1) 
a t 2 a t 2 

where |2#) = -^|0) and |2y) = -/||0) are, correspondingly, biphoton state vectors 
with both photons having either horizontal or vertical polarization, whereas \1h, ly) = 
a^a|/|0) is the state vector of two photons with different polarization, horizontal and 
vertical, a) H and ay are creation operators for photons with horizontal and vertical 
polarizations. The constants 6*1,2,3 obey the normalization condition |C*i| 2 + | C* 2 1 2 + 
I C3 1 2 = 1. In the general case 6*1,2,3 can be complex. But as the phase of the state |^) 
as a whole does not affect any measurable characteristics of qutrits, it can be chosen, e.g., 
in a way making C\ real and positive with the parameters C*2 and C3 being complex (if 
C*i = any other of two remaining parameters C 2 or C3 can be taken real and positive). 
In the general case this leaves four free parameters determining states of qutrits, e.g., 
I C 2 1 and I C3 1 plus their phases <p 2 and </? 3 . 

The main two physically important characteristics of biphoton qutrits are their 
entanglement and polarization. Polarization features of biphoton qutrits have been 
widely investigated since 1999 jU |2] |3] HI |5] . A very interesting and fruitful interpretation 
of these features is related to the description of biphoton qutrits in terms of the 
Stokes vectors on the Poincare sphere (3j 0]. The qutrit state vectors (11. ip were 
considered as reduced to the form of two creation operators acting on the vacuum state 
I 1 !/) = NA'B* |0). Below, such representation of biphoton qutrits will be referred to as 
the "operator factorization", and the operators A' and B\ as "factorizing operators". 
In (3j S] the factorizing operators were considered as functions of angles determining the 
orientations of the corresponding Stokes vectors on the Poincare sphere. The degree of 
polarization of qutrits was found to be related to the angle between these Stokes vectors. 
In particular, maximally polarized biphoton qutris were shown to be represented by two 
parallel Stokes vectors and unpolarized biphoton qutrits, by a pair of counter-directed 
Stokes vectors on the Poincare sphere. 

As to the entanglement of biphoton qutrits, maybe strange enough, but conclusions 
about its existence and degree (if it exists) were often and still remain controversial 
(for comparison, see, e.g., Refs. [6] and [7]). The controversy arises from definitions, 
formulated clearly or not. The simplest definition for pure biphoton states is well 
known. For biphoton qutrits this definition is formulated in terms of the wave function 
^(o~i,o- 2 ) = (cti, cr 2 1 \I>) depending on two polarization variables of two indistinguishable 
photons <7i and o 2 . Each of these two variables can take independently one of two 
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values o~i t 2 = HorV. The state is entangled if the wave function cannot be factorized, 
^(01,02) ^ (p(ai)x{c2) ■ m the frame of this definition and with the obligatory 
symmetry of biphoton wave functions taken into account, the Wootters concurrence 
C [8] characterizing the degree of entanglement of biphoton qutrits was found in recent 
papers P, [10] to be given by 

C=\2C 1 C 3 -C*\. (1.2) 

Moreover, it was found that there is a simple and direct relation between the concurrence 
C and the degree of polarization P for biphoton qutrits, 

C 2 + P 2 = l. (1.3) 

This relation shows that maximally entangled states (C = 1) are unpolarized and, 
oppositely, maximally polarized states (P = 1) are disentangled (C = 0) and their wave 
function is factorized. In particular, as follows from Eqs. ( 11.21) and ( 11.31) . the state 
I If/-, ly) (Ci = C3 = 0, C<i = 1) is unpolarized and maximally entangled (P = and 
C — 1). These and other related conclusions are strongly supported by the analysis 
based on the Schmidt-mode decomposition [91 [TO]. 

In this work we continue investigating the links between entanglement and 
polarization of biphoton qutrits, between the approaches based on the Schmidt 
decomposition and on the analysis of polarization Stokes vectors on the Poincare sphere. 
After a brief reminder about the definitions of the Schmidt modes as the eigenfunctions 
of the reduced density matrix in the next section, in section 3 we suggest a very simple 
procedure of the operator factorization. It shows, in particular, that the operator 
factorization is always doable, and for any given configuration of the qutrit parameters 
C1.2.3 there is only one pair of factorizing operators. In section 4 we describe an 
alternative procedure of the operator factorization based on the use of the Schmidt- 
mode representation. By combining these two approaches we show that there is a way 
of finding the Schmidt modes of a biphoton without solving the eigenvalue-eigenfunction 
problem (section 5). Finally, in section 6 we describe representations of biphoton 
qutrits in terms of the Stokes vectors of either the Schmidt modes or the factorizing 
operators. In this way we establish the full correspondence between approaches based 
on the Schmidt-mode and the Stokes-vector analyzes. 



2. The Schmidt modes 



The Schmidt modes are known to be defined as the eigenfunctions of bipartite 
reduced density matrices. For biphoton qutrits, the derivation is rather simple and 
straightforward [9j[l0]. As mentioned above, the first step is finding the wave function 
^(01,02) corresponding to the state vector (11.11) . Not reproducing here ^ (0-1,0-2) 
completely, let us show only its unfactorable part corresponding to the state |1#, ly): 



^hv(o-i,o- 2 ) = (01, 02 1 Iff, ly) = 



<8> 



• (2.1) 
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The two alternative expressions in the second line of this equation correspond to the 
equivalent functional and matrix forms of writing wave functions. They both contain 
two terms because the wave function of two photons (bosons) has to be symmetric with 
respect to the transposition of variables 0\ # er 2 . In the matrix form the upper and 
lower lines of the columns correspond to the horizontal and vertical polarizations with 
subscripts 1 and 2 indicating photon variables, e.g., as = 5 aij H etc. 

The next step consists of constructing the density matrix p — \& (g) W (Eq. (4.2) of 
Ref. [9]) and its reduction, i.e., taking its trace with respect to one of the variables, e.g. 
variable 2. This yields [9] 



Pr 



( 1^12 IC2I 2 CiCl + c 2 ci \ 
1 ^ 

C1C2 + C2C3 \Q,\2 1 l^ 2 l 



(2.2) 



V y/2 ' 2 J 

As photons are indistinguishable particles and therefore the qutrit wave function has to 
be symmetric, Eq. ( 12. 2 p is also valid for the second reduced density matrix, which can be 
obtained by taking the trace over variable 1. The Schmidt modes are the eigenf unctions 
of the joint reduced density matrix (12. 2p . 

p r ip± = \±(p±, (2.3) 

where A± are the eigenvalues of p r related to the qutrit concurrence C ( II. 2p by equations 

A± = ^(l ± VI - C 2 ) = and C = 2^A+A_. (2.4) 

The Schmidt decomposition for the qutrit wave function has the form 

<P((7i, cr 2 ) = (fi+(o- 1 )ip + (a 2 ) + y/Al y_(«7 1 )y>_(«7 3 ). (2.5) 

Because of the indistinguishability of the photons, pairs of Schmidt modes (</?+, <p + ) and 
(<P-, v 9 -) are formed by identical rather than different functions. This means in fact that 
each of the two terms in the Schmidt decomposition (12. 5p describes a state with two 
photons in the same mode, (+) or (— ) [TT]. For finding the Schmidt modes explicitly, 
one has to solve Eq. ( I2.3p . In the general case such solutions are rather cumbersome. 
An alternative approach is described below. 

The polarization Stokes vector and the degree of polarization can be defined for the 
biphoton qutrit in terms of the reduced density matrix p r as 

S = Tr(p r a) and P = \S\, (2.6) 

where now a is the vector of the Pauli matrices. From here we can deduce immediately 
that the reduced density matrix of biphoton qutrits ( 12.21) coincides with the well-known 
polarization matrix 



Pr - Ppol - o C , „'C ] C 2 • ( 2 - 7 ) 




The polarization matrix is written here in its standard form [12] although with the 
components of the Stokes vector re-numbered. Provided that in the collinear SPDC the 
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biphoton beam propagates along the z-axis with the horizontal and vertical directions 
denoted by x and y, the Stokes parameter 5*3 = 1 corresponds to the linear polarization 
along the x-axis, S3 = — 1 corresponds to the linear polarization along the y-axis, 
Si = ±1 corresponds to the linear polarization along the lines directed at ±45° to the 
x-axis in the xy-pl&ne, and S2 = ±1, to circular polarizations. 



3. Operator factorization 

According to the idea of Ref. [3], the qutrit state vector ( ll.ip with arbitrary parameters 
Ci.2,3 can be reduced to the form of a product of two new creation operators times the 
vacuum state, 

\V) =NA*B*\Q), (3.1) 

where N is the normalization factor such that (^l 1 !/) = 1. The fact that representation 
( 13. ip always exists and is unique, already mentioned in Ref. [3] , can be very easily proved 
as follows. Consider an operator Q such that |^) = Q\0), 

and a second-order polynomial Q(x) of some variable x, associated with Q, 

Q( x ) = ^L x i + C 2 x+^=. (3.3) 
By finding the roots of the equation Q(x) = 0, 



— C*2 i v Cf — 2C1C3 

^ B = • (X4) 

we represent Q(x) as a product of two linear functions of x, 

C 

Q{x) = — 7 ={x-x A )(x-x B )- (3.5) 

As the operators a} H and ay commute with each other, they can be arbitrarily transposed 
and combined as simple numbers. Hence, we can immediately write down the operator 
analog of the algebraic equality (13. 5p . 

C 

Q = ^= (a-H - x A a v ){a ] H - x B a\). (3.6) 

This equation determines the factorizing operators and S' of Eq. ( 13. ip . though it 
leaves undefined their phases y?o an d —fo- 

tf = e*» a *~ XA( ^ and B* = e -*» a * ~ Xb< ^ , (3.7) 
' + \xa\ 2 \/l + \xb\ 2 



The phase factors in these definitions do not affect usual commutation rules for the 
photon annihilation and creation operators in a given mode, 

[A, A 1 ] = AA* - A^A = 1 and [B, fit] = 1. (3.8) 
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But they do affect the cross commutators [A, fit] and [fi, At] = [A,fit]t. As shown in 
the following section, these commutators are always real, whereas the parameters xa,b 
(13.4p can be complex, xa,b = \ x a,b\ vx${iip a,b) , where <£a and ips are phases of xa and 
xb- Thus, the phase ipo can be found just from the condition that the cross-commutators 
are real 

= e^l + x^) = U+ 

^l + WKi + W) J(i + \x A \*)(i + \x B \i) 



Explicitly, the phase <po is given by 
1 

ifo = - tan 





xa\ 


\x b 


sin(v9 B - <p A ) 


1 + 


xa\ 


\xb 


cos((p B - <Pa) 



(3.10) 



The commutator (I3.9P is a very important characteristics of the biphoton qutrit states 
and, in particular, it determines the normalizing factor N in Eq. ( 13. ip 



iV= . 1 . (3.11) 

4. Operator factorization in the Schmidt representation 

The Schmidt decomposition (12 .5p provides an alternative way for the operator 
factorization of qutrit state vectors. In this approach we can further specify the earlier 
introduced parameters, such as the commutator of the operators A and fit ( 13. 9p . as well 
as the definition (13.71) of the operators At and fit themselves. 

As mentioned above, owing to the symmetry of the qutrit wave functions, the 
Schmidt decomposition ( 12. 5 p presents qutrits as a coherent superposition of two 
orthogonal two-photon states, |2 + ) and |2_). For this reason, directly from Eq. (12 .5^ 
we can find the operator Schmidt decomposition of the qutrit state vector, 

— f«/I7 nX 2 4- ,f\~ o.t : 



m= V2 [v X+a+ +V A - a -JI°>> ( 41 ) 

where a + and a_ are the photon creation operators for the Schmidt modes (+) and (— ). 
In terms of a}± the operator Q (13. 2p takes the form Q = (\/A~+ a) + + \/X~ at ) / \/2, and 
it is easily factorized to give the above-defined and appropriately normalized operators 
A' and fit in the Schmidt representation, 



At 
fit 



Ai /4 al + iX'l'al \Tal + iX]/ 4 al 



1/2 - (1 + C y/A ' 

(4.2) 



xY 4 al_ — iXY^aL AV 4 at — iXY^at 



/ A7+ a/A 



1/2 " (! + c y/4 



The commutator of the operators A and fit can be now written in terms of the 
Schmidt eigenvalues X±, or the concurrence C, or the degree of polarization P, 



[A, fit] = [fi, ^] = ^-^Z = J?—£ = £ , (4.3) 

VK + t/xz Vi + c i + vi - p 2 
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and similarly for the normalization factor in Eq. (13. ip . 



N 



V2 



c 



(4.4) 

Eq. ( 14. 3p shows that in the maximally entangled unpolarized states (C = 1, P = 0) the 
commutator [A, B^] equals zero, which means that the operators A> and B^ represent 
orthogonal modes, as, e.g., a} H and a v . In the case of disentangled maximally polarized 
states (C = 0, P = 1) the commutator [A, £?'] equals 1, which indicates that the 
operators ^ and B^ coincide with each other. Eq. (14. 3[) can be inverted to express the 
concurrence C via the commutator [A, B^], 

In principle, this formula can be considered as a definition of the concurrence alternative 
to that of Eq. ( 11.21) . Though it is hardly easy to find the commutator [A, B'] before 
finding the concurrence, Eq. (14. 5p is important because it shows that for getting 
disentangled states of biphoton qutrits it is not sufficient to provide the operator 
factorization. It is necessary to have the commutator equal to unity, [A, St] = 1. 
In other words, the operator factorization by no means leads automatically to the 
factorization of the biphoton wave function and, hence, the operator factorization does 
not guarantee disentanglement. The simplest and most often discussed example is the 
state \1h, lv)- The operator factorization is present automatically with the factorizing 



operators A 1 



and St 



But the commutator [A, B^\ = [a#, a v ] equals zero 



and hence, there is no factorization in the biphoton wave function ( 12. ip . the state is 
maximally entangled and unpolarized. 



5. Finding the Schmidt modes without solving the eigenvalue-eigenfunction 
equation flUS) 

Eqs. ( 14.2p can be inverted to express the Schmidt-mode creation operators a) + and at 
in terms of the factorizing operators A^ and B^ 

1/4 




(At + St), 



(5-1) 



^ -At). 



By substituting here expressions (13 ,7p for A^ and B^ via a) H and ay, we finally get the 
general expressions for the Schmidt-mode creation operators, 

1/4 



1 fl + C 



e lLpo x A 
1 + \xa\ 



-vp 



i + M 2 a/i + |x b | 2 



-upo 



+ 



XB 



\Xb\ 



(5.2) 
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and 



1/4 



i + M 2 Ji + \x A \ 2 



-upo 



X B 



X A 



1 



\Xb\ 



FA 



(5.3) 



These operators can be checked to obey the usual commutation rules for the creation 
and annihilation operators in two orthogonal modes, 



.,a\] = 1, [a_,aL] = 1, [a+,at] = [a_,a+] = 0. 



(5.4) 



<p±{a) = (a\ai\0) = 



(5.5) 



The Schmidt modes corresponding to the operators a+ (15 .2D and at (I5.3P are given by 

'(i H \ai\oy 

H«±|o) 

It should be stressed that the Schmidt modes are found here with the help of 
simple algebraic operations and without solving the eigenvalue-eigenfunction equation 
(12. 3p . which is rather unexpected. 

As an example, let us consider a simplified biphoton qutrit of the form 

|\&) = N a} H (cosaajj + sin a ay) |0). (5.6) 

This is the case when the qutrit state vector "automatically" has the form NA'B'\0) 
and no factorization procedure is needed. The factorizing operators A 1 and are then 
given by 

A* 



a* H and 



cos aa{f + sin a ay. 



Their commutator and normalization constant iV are equal to 
[A, B ] ] = cos a and N 



cos 2 a 



In the original form (II. lft . the parameters of the qutrit ( 15. 6ft are equal to 

V2 cos a „ sin a 



C x 



Co 



cos 2 a ' 



C, = 0. 



(5.7) 
(5.8) 

(5.9) 



\/T + cos 2 a' a/1 + 

These parameters correspond to the following concurrence C, degree of polarization P, 
and the eigenvalues of the reduced density matrix A±: 

2C0Stt (5.10) 



C 



sin 2 a 



cos 2 a ' 



X- 



1 



1 ± VT^c< 



P = VT 
1±P 



C 2 



cos 2 a 



1 (1 ± cos a) 2 

2 1 + cos 2 a 



(5.11) 



As the factorizing operators (15. 7p of the qutrit (I5.6P are known "automatically", the 
Schmidt-mode creation operators can be found directly from Eqs. (15. ip rather than 
from more complicated equations (15. 21) and ( 15.31) . 

1 



2(1 + cos a) 

i 

2(1 — cos a) 



cos a a 



sin a a 



- 1 + cos ct) oJjj + sin ct a v 



(5.12) 
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From these equations and with the help of Eqs. (15.51) we easily find explicit expressions 
for the Schmidt modes of the biphoton qutrit (15. 6p 

= 1 ( 1 + cos a \ _ / cos(a/2) \ , 

'2(1 + cos a) \ sina ) \ sm(a/2) ) ' 

V- = , % ( ~ l + C ° Sa 1 = i ( ~ Ma l^ ) • (5.14) 

2(l-cosa)V sma / V cos(a/2) 



6. Poincare-sphere representation 

6.1. The Stokes vectors of single-photon states 

Following earlier works [3J H] let us consider representation of the above-derived results 
in the space of Stokes vectors on the Poincare-sphere. Let us start with the general 
expressions (13.71) for the factorizing operators and B'. These operators generate a 
pair of single-photon states and the corresponding wave functions: 

|1 A ) = At|0) = . 1 (4, - x A al)\0), (6.1) 
|l B ) = fft| )= . 1 {a} H - x B a} v )\Q) (6.2) 



x B 



and 



= I = ■ ( 1 ),1/>B= 1 = : ( 1 ) , (6.3) 

where, as shown above, the parameters xa.b can be expressed via the original qutrit 
parameters by equations (13.41) . We have dropped the phase factors occurring in operators 
A' and 5' (3.7) because they do not affect the polarization matrices of the states 
considered below. 

The density matrices of the single-photon states \1 A ) and \1 B ) are easily found from 
their wave functions to be given by 

pA = TXT — 12 f 1 7%)> ( 6 - 4 ) 

and the same for p B with the substitution xa — > x B . By identifying the density 
matrices of the states \1a) and \1 B ) with the corresponding polarization matrices p^j 
and pj^j (having the same form as the polarization matrix in Eq. (12. 7p ) we can find the 
components of the Stokes vectors and of the states \1 A ) and \1 B ): 

Si A) = (l-\x A \ 3 )/(l + \x A \*)=<x*0 At 
S[ A) = -2Re(x A )/(l + M 2 ) = sin^ A cos^> A , ( 6 - 5 ) 
= — 2lm(x A )/(l + \xa\ 2 ) = sin6U sinip A , 
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S [ 3 B} = (l-\x B \ 2 )/(l + \x B \ 2 ) = cos8 B , 
S{ B) = -2Re(x B )/(l + \x B \ 2 ) = sm8 B cos<p B , (6.6) 
S2 = — < 2lm(x B )/(l + \x B \ 2 ) = sin8 B simp B , 

where 8 a and 8 B are the angles between the vectors and and the x-axis on the 
Poincare sphere (FigHJ). The values (fA and (f B are the angles between the horizontal 
axis 045° _L Ox on the Poincare sphere and projections of the Stokes vectors and 
S ^ on the plane perpendicular to Ox. 




Figure 1. The Stokes vectors (with spherical coordinates 8a, ^a) and 

(|6.5I6.6J) . as well as the biphoton Stokes vector S Uph (j6TT|) . (pl9| . (jO0|) . on the 
Poincare sphere. 



As it follows from Eqs. (16 .5D and (16. 6\\ . the absolute values of the Stokes vectors 



and are equal to unity, = 
8ab is determined by their scalar product, 

cos 8 AB = S {A) ■ S iB) = cos 8a cos 8 b 



1. The angle between these vectors 



sin 8a sin cos^ipA — Vb)- (6.7) 



Schmidt-mode representation and Poincare- sphere images of biphoton qutrits 



Let us analyze now the Poincare-sphere image of the general- form biphton qutrit (13. ip . 
One possibility of its investigation consists of finding from this equation the original 
qutrit's parameters Ci i2 ,3, substituting them into the general expression for the reduced 
density matrix (12. 7ft . and finding from this last equation the biphoton Stokes parameters, 
which involves rather cumbersome transformations and formulas. An alternative, more 
elegant and even more informative approach is based on the use of the Schmidt-mode 
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representation. In this representation the Schmidt modes and their single-photon density 
matrices are very simple: 

The Stokes vectors 5+ and 5_ corresponding to the Schmidt modes ( 16.81) have 
components 

4+) = 1, S { 2 +) = S{ +) = and St ] = -1, 4 _) = = 0. (6.9) 

In terms of the Poincare sphere this means that the transition to the Schmidt-mode 
representation is realized by means of a rotation of the Poincare sphere that makes 
the Stokes vectors of the Schmidt modes directed along the former (H, V) axis. Then 
this axis turns into the (+, — ) axis, so that the Stokes vectors and point, 
correspondingly, at the positive and negative directions along this axis (Fig. E}. 
The biphoton reduced density matrix (12. 7p is also very simple in the Schmidt-mode 




Figure 2. a: Poincare sphere and the Stokes vectors in the Schmidt- mode 
representation; b: A scheme of two Schmidt-mode two-photon states 



representation 

/ \ . n \ 

(6.10) 




As well as the one-photon Schmidt-mode Stokes vectors S^^, the biphoton Stokes vector 
gbiph on ly one non-zero component, 

S ^ h = \ + -\_ = p W ith S? ph = S? iph = 0. (6.11) 

On the Poincare sphere the vector S htph is directed along the same axis (+, — ) as the 
Schmidt-mode Stokes vectors S^> . As usual, the absolute value of S b%ph coincides with 
the degree of polarization per photon P for biphoton polarization qutrits [31 H] , 



S^ ph = P. (6.12) 
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Finally, the factorizing operators and are given by superpositions of the 
Schmidt-mode creation operators (14. 2J) . In the Schmidt-mode representation the one- 
photon states generated by these operators, A^|0) and-B^O), have the following wave 
functions 



A 



1/4 



(1 + C) 1 / 4 I iX ]i 



1/4 



A 



1/4 



1/4 



and density matrices 



PA.B 




Components of the Stokes vectors of these states are S[ A ^ 



1/4 



and 



si B) 



VT+C 
2(A + A_) 1 / 4 




(6.13) 



(6.14) 



(6.15) 



VT+C 

Thus, in the Schmidt-mode representation the Stokes vectors and are located 
in the vertical plane containing the horizontal axis (+,—), and their directions are 
symmetric with respect to this axis. The angles between 5 ,( - j4 ' B ^ and the biphoton Stokes 
vector S tnph \\(+, —) are determined by S^" 4 -projections of the vectors S^ A ' BS> on the 



cos (8 AB /2) 



l-C 



p 



[A,B\ (6.16) 



V 1 + C 1 + yfl^Pl 

where Qab is the angle between the vectors and S^ B \ 

By inverting the relation between cos (#ab/2) arid the concurrence C, we can express 
the latter in terms of the angle Qab between the vectors and S^: 

1 — cos 2 (#ab/2) 1 — cos 8 ab 



C 

1 + cos 2 (#ab/2) 3 + cosW 
And the expression for the degree of polarization P via 9ab has the form 

2a/2(T 



p = vi-c 2 



cos 6 'ab) 4|cos(0ab/2)| 



(6.17) 



(6.18) 



3 + COS Q AB 3 + COS Q AB 

This result for the degree of polarization agrees perfectly with that of the paper [I] . 

As the biphoton Stokes vector is directed along the bisector of the angle between 
the Stokes vectors and S^ B \ in the vectorial form the biphoton Stokes vector can 
be written as 



5 



biph 



s 



biph 



§(A) + g(B) 



2cos(^ B /2) 
(3 + cos 6ab)(1 + cos Qab) 



P- 



+ cos Bab) 



2(1 



(6.19) 
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Figure W? illustrates an alternative representation of a qutrit in terms of the 
Schmidt-mode Stokes vectors and S^'. The picture shows pairs of these Stokes 
vectors, oppositely directed and having the weight coefficients A + and A_. The total 
biphoton Stokes vector equals the sum of all four Stokes vectors of the Schmidt modes 
with the weight coefficients, and half of this sum gives biphoton Stokes vector per one 
photon: 

gbiph = A+ £(+) + A _£(-) _ ( 6>20 ) 

The results of this section are invariant with respect to rotations of the Poincare 
sphere. The only changes concern the orientation of the (+) — (— ) axis and the plane 
containing the Stokes vectors Sa and S B . But at any orientation of the (+) — (— ) axis 
and (Sa, Sb) -plane, the Stokes vectors of the Schmidt modes S^\ of the factorizing 
operators S^ A ' B \ and of the biphoton state as a whole S btph belong to the same plane. 
The biphoton Stokes vector is directed along the same direction as S^\ and this 
direction coincides with the bisector of the angle 9ab between Sa and S B . Eqs. ( 16. 1 7ft - 
( I6.20p remain valid without any changes in any frames with rotated Poincare sphere. As 
to other derived equations (16. 9p . (16. lip . (I6.12p . (I6.15p . and ( I6.16p . they remain valid too, 
though the components of Stokes vectors S3 and S2 have to be understood as projections 
on the directions of the turned (+) — (— ) axis and turned \Sa, Ss^-plane, rather than, 
e.g., projections on the V—H and R—L axes of the standard Poincare-sphere orientation 

of Fig. m 

7. Possible experiments 

A possible experiment on the selection of polarization Schmidt modes of a qutrit is 
shown in Fig. [3^l. The qutrit state is sent to a polarization beam splitter (PBS) preceded 
by a quarter-wave and a half-wave plates. The plates are oriented in such a way that 
photons in the polarization Schmidt mode \^+) become horizontally polarized p3]. Then 
the orthogonally polarized photons in the \<p~) mode automatically become vertically 
polarized. The qutrit state vector after the plates is transformed into a weighted 
superposition of photon pairs in vertical and horizontal modes, 

\*) = y /^\2 H ) + e H +fil\2 v ), (7.1) 

where <fi is the relative phase of the Schmidt modes |</?_) and \(p+) |llj . 

After the PBS, the pair |2#) goes into the transmitted output port and the 
pair |2y), into the reflected output port. Each photon pair can be detected as 
a coincidence of single-photon detector 'clicks'. For registering such coincidences, 
beamsplitters (BS) followed by pairs of detectors are introduced in both output ports 
of the PBS. Then, coincidences will be observed between the counts of either detectors 
Di + , D2+ or detectors -Di_, -D2-, and their rates R± will scale as the Schmidt eigenvalues 
X± = R C ± /(R C + + R C J). But no coincidences between the outputs of the detectors in 
different ports of the PBS will be observed. The polarization Schmidt modes can be 
found experimentally by finding positions of the quarter-wave and half-wave plates at 
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which the coincidence counting rates between the detectors in different output ports of 
the PBS turn to zero. This approach is most appropriate if the qutrit parameters Ci 2 3 
are not known in advance. Further, from the experimentally found parameters A±, with 
the help of Eqs. (I2.4p one can easily calculate both the qutrit concurrence and degree 
of polarization. 

In the absence of detectors' dark counts and parasite light, the beamsplitters in the 
PBS output ports are not necessary and only a single detector can be used in each port. 
Note that this method of finding A± can be made invariant to the quantum efficiencies 
of the detectors. Indeed, instead of comparing the counting (coincidence) rates at the 
two output ports of the PBS one can compare the counting (coincidence) rates in the 
same port but at two different positions of the plates, corresponding to the Schmidt 
mode \ip+) or sent to the same port. 

It is worth mentioning that the four parameters determining a qutrit can be also 
introduced as the orientations of the plates transforming the qutrit from the form (14. ip 
into the one (17.11) . one of the Schmidt eigenvalues, and the phase <p. Experimentally, the 
latter can be measured in the same scheme as shown in Fig. [3^ but with the PBS turned 
for 45°. Indeed, if the number of qutrits at the input of the PBS is N, and the quantum 
efficiencies of detectors D x+ and D 2+ are, respectively, 771 and r/ 2 , then the numbers 
of their coincidences with the PBS oriented at 0° and 90° will be Rq = r)iT]2\+N/2 
and -R90 = f]ir]2\-N/2, respectively. At the same time, for a 45° orientation of the 
PBS the number of coincidences will be i? 45 = 771772(1 + 2 y/\+ A_ cos 2<p)N/ 4. From the 
measurement of these three parameters, one can infer both the Schmidt eigenvalues and 
the phase <f) without knowing the quantum efficiencies and the initial number of pairs: 

, Ro , R90 n , 2i?4 5 — Rq — Rg , . 

A + = — — - , A_ = — — — , cos 2<p = , . (7.2) 

Ro + R90 Ro + R90 2V-R0-R90 




Figure 3. a: Experimental setup for observing polarization Schmidt modes of a 
qutrit \^>); b: Experimental setup for demonstrating polarization entanglement of a 
qutrit I*). 

For the physical interpretation of the polarization entanglement of a qutrit, one can 
suggest another experimental scheme (Fig. |3b). The qutrit after the polarization 
transformation, is directed to a 50% non-polarizing beam splitter (NPBS). Half of the 
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photon pairs will exit into a single output port of the NPBS, but the other half will be 
split between the two ports. This part of the state can be written as 

|tf> 12 = Ih2) + e 2i<t> y/\l \1 V1 , l y2 ), (7.3) 

where the indices HI (H2) denote horizontal polarization in the transmitted (reflected) 
beam, and similarly for the vertical polarization. 

This state has been considered in Ref. [H], and it was shown to be maximally 
entangled only if the Schmidt eigenvalues are equal, A + = A_ . This exactly corresponds 
to the case where C = 1 and the qutrit |^) is fully entangled. If A + ^ A_, the state is 
non-maximally entangled, and its degree of entanglement can be measured, as shown in 
Ref. [Hj . from the orientations of the polarizers P\, P 2 at which the coincidence counting 
rate for detectors D\, -D 2 vanishes. 

8. Conclusion 

We have considered the general case of a biphoton qutrit, physically represented by 
a photon pair in a single frequency and wavevector mode. We have found a simple 
procedure of factorizing the qutrit state vectors, i.e., representing them as products of 
two single-photon creation operators in certain polarization modes acting on a vacuum. 
The procedure clearly proves the existence and uniqueness of such a representation. 
We have shown that the described operator factorization of state vectors is not related 
to the factorization of the biphoton polarization wave function and does not exclude 
a possibility of entanglement for biphoton qutrits. The degrees of their entanglement 
and polarization are found to be related directly to the commutation features of the 
operators factorizing the state vectors by Eqs. (14. 3ft and ( 14.51) . Further, we considered 
the polarization Schmidt decomposition of a generic biphoton qutrit and derived it 
both using the standard method and the factorization procedure. For the qutrit and 
its Schmidt modes, the Stokes vectors have been calculated and shown on the Poincare 
sphere. Finally, an experiment has been suggested to find the polarization Schmidt 
modes of a qutrit and to measure its degree of entanglement and polarization. A physical 
interpretation is proposed for the polarization entanglement of a biphoton qutrit. 
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